Amplifying single impurities immersed in a gas of ultracold atoms 
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We present a method for amplifying a single or scattered impurities immersed in a background gas 
of ultracold atoms so that they can be optically imaged and spatially resolved. Our approach relies 
on a Raman transfer between two stable atomic hyperfine states that is conditioned on the presence 
of an impurity atom. The amplification is based on the strong interaction among atoms excited to 
Rydberg states. We perform a detailed analytical study of the performance of the proposed scheme 
with particular emphasis on the influence of inevitable many-body effects. 



Impurities inside interacting quantum fluids provide 
a variety of interesting physical effects in fields such as 
high-T c superconductivity or plasma physics. Recent 
experimental advances in creating and controlling im- 
purities in ultracold quantum gases have demonstrated 
that these systems are ideally suited for studying impu- 
rity phenomena such as transport in strongly correlated 
lD-Bose gases pQ, polaron physics in highly imbalanced 
Fermi-gas mixtures [2], or Rydberg-excitation saturation 
in Bose- Einstein condensates [3] . Since detection of single 
(neutral) impurities in a quantum gas poses severe exper- 
imental challenges, these experiments all rely on measur- 
ing either ensemble averages or effects of the impurities 
on the bulk host medium. Here, we develop a novel single 
shot detection scheme for single impurity atoms, based on 
strong amplification of the detectable signal by convert- 
ing the impurity into an internal state change of a large 
number of nearby background atoms as shown in Fig. [T^,. 
Our scheme is applicable to any form of neutral impurity 
in an ultracold gas and can be implemented with stan- 
dard absorption-imaging techniques available in virtually 
every existing experiment, e.g. one possible application 
is the visualization of Rydberg crystals in a Bose gas [3] . 

In the experimental situation we have in mind the 
background gas is initially polarized with all atoms be- 
ing in a stable hyperfine state \A). The average distance 
i?b between the impurities, which can be atoms of the 
same species in other internal states or a different atomic 
species altogether, is much larger than the average dis- 
tance between background atoms. As illustrated in Fig. 
[TJd, our amplification scheme performs a conversion of 
background atoms from \A) to a second hyperfine state 
\B) only within a sphere of radius R' c centered at the im- 
purity. These atoms can be detected via state selective 
absorption imaging in a single shot and provided that 
i?b > 2R' C a discrimination of the positions of the im- 
purities can be achieved. We show that our scheme can 
be directly implemented in current experiments with ul- 
tracold Rubidium gases yielding amplification factors of 
~ 10 — 100 atoms per impurity. 

We start from a setting where all impurity atoms have 
been electronically excited to a Rydberg n'S-state. Such 



FIG. 1. a: After the amplification step the atoms of the back- 
ground gas have changed their internal state [from \A) (green) 
to \B) (orange)] in the vicinity of an impurity (black). State 
selective absorption imaging of the background gas then per- 
mits the spatial discrimination of the impurity atoms, b: 
Length scales. The mean distance Rb between the impurity 
atoms is much larger than the mean interparticle distance in 
the uniform background gas. Atoms within a sphere with 
radius R' c change their internal state. Background atoms lo- 
cated within R c are interacting (see text for explanation). 

situation can be easily achieved by a resonant 7r-pulse 
from the ground state. As an example, one can con- 
sider the case of strongly interacting Rydberg atoms in 
a cloud of ultracold Rubidium atoms, where the typical 
distance between impurities is given by the so-called Ry- 
dberg blockade radius [5J [B] which can be on the order 
of ten micrometers [Tllllj. To effectuate the conditional 
transfer between the hyperfine states of each background 
atom we use the atomic levels depicted in Fig. [2^,. This 
scheme has been thoroughly discussed in Ref. ]T2] in the 
context of a mesoscopic quantum gate. The level \P) is 
a metastable low-lying P-state over which the hyperfine 
states of interest (\A) and \B}) are off-resonantly cou- 
pled (detuning A p , Rabi frequency f2 p ). The state \P) 
is furthermore coupled to a Rydberg nS state \R) with a 
laser of Rabi frequency fi c whose frequency is chosen such 
that the states \A) and \R) are in two-photon resonance. 
The interaction between the Rydberg state and the im- 
purity, denoted by V, acts as an effective detuning for the 
\R) state. Thus the Hamiltonian of a single background 
atom is given by H schcmc = fi p [\A) (P\ + \B) (P\ + h.c] + 
n c [\R)(P\+h.c] ~ A P \P)(P\ + V\R)(R\, with h = 
1. Upon adiabatic elimination of the far-detuned P- 
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state (|A p | 3> Sip, £l c ) and the introduction of the su- 
perposition states |t> = (l/Vf)[\A) + \B)] and ||) = 
(l/v2) [\A) — \B)] the system reduces to a two-level sys- 
tem plus an uncoupled third level (see FigJS]}) . The cor- 
responding Hamiltonian reads 

H ad = [|t) (R\ + h-c] + (A + V)\R) (R\ - e \i) (i\ 

with the (effective) Rabi frequency and detuning il — 
\/2n p Q c / A p and A = (fi^ — 2f2p)/A p , respectively and 
where e = 2f2 p /A p is the energy of the |f)-state. This 
scheme formally corresponds to an off-resonant Raman 
coupling between \A) and \B) via the Rydberg state. 
Note however, that here 51, A and e are not independent 
from one another. 

To illustrate the mechanism behind the amplifier, we 
consider a single impurity atom localized at the origin 
of the coordinate system. The interaction energy be- 
tween this impurity (excited to the Rydberg state n'S) 
and the fc-th atom of the background gas located at a 
distance |R/-| = Rk is given by 14 = Cg/i?|. Here 
Cg is the van-der-Waals coefficient for the n'S — nS 
interaction. This interaction energy together with the 
detuning A defines a spatially varying effective detun- 
ing of the Rydberg state of the fc-th background atom: 
A fc = A + V k = A(l + R' c 6 /R 6 k ) where we have defined the 
radius R' c = (Cg/A) 1 / 6 . As we will show, this spatially 
dependent energy shift is the origin of the conditional 
transfer between the hyperfine states of the background 
atoms. For £! c ^> Sl p , which we assume throughout, the 
probability to excite a Rydberg state is ~ (Sl/A) 2 <C 1. 
We can thus adiabatically eliminate the state \R) and 
find that each background atom evolves effectively under 

the Hamiltonian H n = -e fl/6 fl | fi g |t)(tl + 11X11 • We 
switch on the effective laser coupling (Rabi frequency $1) 
for a time = 7rA p /2f2p = 7r/e (7r-pulse). This results 
in a conditional transfer between the hyperfine states as 
a consequence of the spatially dependent phase difference 
that is accumulated between the states |t) and \\) under 
the action of H^: If the background atom is close to the 
impurity (Rk — > 0), the coupling realizes the transition 
\A) —> \B), i.e. the atoms of the background gas change 
their state since |t) acquires a phase of ir relative to |f). 
Far from the impurity (Rk 3> R' c ) the acquired phase 
difference is zero and hence all atoms remain in the state 

While the mechanism is intuitive for a single back- 
ground atom, its generalization to the many-body case 
is not straightforward. This is due to the fact that the 
atoms forming the background gas are coupled to Ry- 
dberg states during the amplification procedure which 
then leads to an interaction among themselves. The im- 
portance of this interaction has been impressively demon- 
strated in recent experiments [T3] and theoretical work 
[TH - TTT] that investigate electromagnetic transparency in 
Rydberg gases. In the following we give a detailed ac- 
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FIG. 2. Internal level scheme of the atoms of the background 
gas. a: The two hyperfine ground states \A) and \B) are 
coupled off-resonantly to the intermediate \P) state by means 
of a laser (detuning A p , Rabi frequency fi p ). A strong laser 
with Rabi frequency f2 c couples \P) to a Rydberg state \R). 
The interaction between the impurity and the background 
atom V gives rise to an effective detuning of the Rydberg level, 
b: The scheme can be mapped into a two-level system formed 
by It) = (1/V2) [\A) + \B}] (with energy e) and the Rydberg 
level, coupled with an effective Rabi frequency SI and detuning 
A + V and an uncoupled level j|) = (1/V2) [\A) - \B}]. 



count of these many-body effects. The interaction energy 
between two background atoms in the Rydberg state is 
given by V km = C 6 /|R fc - R m | 6 , where R fe and R m are 
their respective spatial positions and C§ is the van-der- 
Waals coefficient for the nS — nS interaction. We proceed 
by adiabatically eliminating all many-body states that 
contain one or more Rydberg excitations [HI US] . After 
this elimination, we can rewrite the Hamiltonian in terms 

(k) 

of the spin operator Sz that acts on the background 
atom located at R^ and we are left with the effective 
Hamiltonian 

H = J2 ^4 fe) + E MkmSW + M km S^S^ (1) 
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This Hamiltonian is valid for atomic densities p that obey 
p(C 6 /A) D / 6 (n/A) 2 < 1, where the two-body potential 
dominates the many-body interactions |20) . with D be- 
ing the dimension of the system. The Hamiltonian ([!]) is 
interpreted as follows: The first term is the sum of in- 
dependent atoms that undergo a rotation about the spin 
z-axis depending on the distance to the impurity (similar 
to the Hamiltonian i/j-4.). The second one effectuates an 
extra rotation about the z-axis due to the interaction be- 
tween the Rydberg states of the background atoms. The 
third term is a spin-spin interaction. We first focus our 
discussion on the first two terms, and show later that the 
third term can be understood as a squeezing interaction. 

We employ a continuum representation for the Hamil- 
tonian such that J^k ~* P I ^ D R- The double sum 
in the second term then turns into a double integral 
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over the spatial coordinates Ri and R2. To simplify 
the Hamiltonian we perform a coarse graining procedure 
in which we describe all atoms located at distance R 
from the impurity in terms of a spin density S Z (R), i.e. 
Y, Rk =R S z = dRR D - 1 S z {R) JdfiD-i EQ- This is mo- 
tivated by the fact that in a background gas with uni- 
form density the interaction energy of each atom of the 
gas can only depend on its distance R to the impurity. 
We introduce the relative and center of mass coordinates 
r = Ri — R 2 and R = (Ri + R2V2, respectively, inte- 
grate over the relative coordinates and average over the 
solid angle of the center of mass. This leads us to 

H= (dRR D - l S z {R) J dtl D - 1 [h{D,R) + f 2 (D,R)] 



with 
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where x k = R' c 6 /R 6 k , x = R' c 6 /\r\ 6 and a = C 6 /C' e . One 
can rewrite the previous expression as 



A^ 



f 2 (D,R) = ^ D pR?g D (R,a), 



(2) 



where we have defined the radius R c = (Cq/A) 1 / 6 and 
where £d is a numerical factor that depends on the di- 



mension D of the system: £1 



and 



32 l/6 ! ^ — 3 ^32 1 /3 

£3 = 3^72- The function gu(R,a) is for R — > and 
1 for R ^> R' c . It only weakly depends on the dimen- 
sion and a, and it is well approximated by gn(a,R) ~ 
R 7 /(R 7 + lOR'J), provided that < a < 10. 

Next we discuss the actual dynamics of the amplifica- 
tion procedure in the many-body picture. Initially the 
background gas is prepared in the state 1^4)^, where 
N = f dn D ^xJ dRR D ~ 1 p. The time-evolution of this 
state can be studied for each shell of radius R indepen- 
dently. The mean density of |B)-atoms in such a shell is 
given by 



p B (R,t) = ^{l 



cos[t(f 1 (D,R) + MD,R))]}. (3) 



We are interested in the state of the background atoms 
after a 7r-pulse of duration r w ~ n/e. Close to the impu- 
rity (R ->■ 0) we find Pb(0, t„) — p while for R^> R' c the 
density is approximately 
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FIG. 3. a: Background density of \B) atoms after the Raman 
transfer versus the parameter pR c . The data shown corre- 
spond to numerical simulations of a one-dimensional lattice 
gas of N = 20 atoms. For each value of £l p /£l c = 0.1 (di- 
amonds) 0.15 (squares) and 0.2 (circles), two cases are de- 
picted: the data using the complete Hamiltonian |l| (red 
dotted lines) and without taking into account the third term 
proportional to SzSz (blue solid lines), b: Density of 
atoms in the \B) state for the same system, with f2 p /f2 c = 0.1 
and one impurity atom at each end of the lattice (blue dots). 
The values of the parameters considered are: principal quan- 
tum numbers of the background and impurity atoms n = 39 
and n' = 62, respectively; n p /(27r) = 10 MHz; Q c /(2ir) = 
100 MHz; A p /(2tt) = 1GHz; p = 2^m _1 . The impurity is 
amplified by a factor ~ 10. 



where we have expanded in f 2 (D, R)t^ <C 1. This back- 
ground of |_B)-atoms is caused by a mean field energy 
shift due to virtually excited background atoms that in- 
troduces a small phase difference between the states |f) 
and 1 4-). This clearly illustrates the importance of the in- 
teraction between the Rydberg states of the background 
atoms. This effect can be minimized by choosing either 
Rydberg states in which the background atoms interact 
weakly (small R c ), by decreasing the ratio Q p /Q c or by 
reducing the density p. The latter is not desirable as this 
reduces the amplification factor, i.e. the total number of 
I B) -atoms produced by an impurity which is given by 



J dil D -! J dRR?- 1 \p B (R,T v )-pB(pO,T lr ) 

P R' c D Jdn D ^ 



tt 2 

7Z3 + X D TT-^^DpRc 



(5) 



where the second line is again the result of an ex- 
pansion in f2(D,R)r 7r <C 1. The parameters -fo = 
{1.02,0.53,0.39} and \ D = {0.16,0.24,0.41} have been 
calculated with the analytic approximation to go(o:,R). 

Let us discuss the role of the last term (spin-spin inter- 
action) of the Hamiltonian ([!]) which we have omitted so 
far. Its effect can be estimated as follows: Each atom of 
the background gas is interacting with other background 
atoms in the volume ~ £dRq 9d{R, ct). Within this vol- 
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FIG. 4. Ratio a = Ce/Cg for pairs of Rydberg S-states of 
Rubidium with principal quantum numbers (n, n') that ac- 
complish < a < 1. In the inset, the value of the C§ in 
atomic units calculated by using second order perturbation 
theory as in Ref. [24] . 



ume (see Fig. [TJ)) only one of the background atoms 
can be (virtually) excited. We consider now a picture in 
which the space is divided into independent regions of the 
latter volume. This allows us to rewrite the last term of 
Hamiltonian ([I]) as a spin squeezing interaction |22j that 
in general affects the transfer of background atoms from 
\A) to \B). This effect can be approximately accounted 
for by multiplying the factor co^ DpR ? 9o{R,a)-l (^tj in 

front of the cos-term in eq. ([3| (see Ref. [23]). Due to 
the functional form of go {R, a) , this factor is close to one 
near the impurity and atoms are transferred between the 
hyperfine states as anticipated. The effect for R ^> R' c 
is demonstrated in Fig. [3^, which shows results obtained 
from numerically solving Hamiltonian ([lj on a ID back- 
ground lattice gas of 20 Rubidium atoms. As predicted 
the background density of |_B)-atoms pb(oo,t w ) grows 
with increasing pR c . The effect of the spin-spin inter- 
action is marginal. Fig. [3}} shows that two impurities 
at a distance i?b ~ 10 pm can be resolved by optically 
imaging atoms in the state \B). The agreement with the 
analytical prediction is excellent. 

We now present experimental parameters for the im- 
plementation of our scheme in current experiments in ul- 
tracold gases of Rubidium. In order to maximize the con- 
trast, we need the background atoms to interact strongly 
with the impurity atoms and weakly among themselves, 
i.e. a = Cq/Cq < 1. Moreover, the sign of the two in- 
teractions must be positive in order to avoid unwanted 
resonances. This is accomplished by only a few combi- 
nations of principal quantum numbers (n, n') as shown 
m Fig. [4j To avoid loss of atoms from the intermediate 
\P) state, the duration tv of the Raman transfer must be 
shorter than the lifetime of the state (r p pa 27 ns for Ru- 
bidium 5P state) divided by the probability of exciting 
the state (given by f2p/Ap). This yields a lower limit for 
the detuning A p /(27r) > 50 MHz. Moreover, 7V has to (2010) 



be shorter than the lifetime of the Rydberg state of the 
impurity, typically on the order Tr pa 100 /is. This sets 
a minimum value for ^ p /(2tt) > ^7r(A p /(27r))/(2T R ), 
that will depend on the particular choice of the impu- 
rity state. Finally, £l p /fl c <C 1 for the Hamiltonian 
([!]) to be valid. Taking all this into account the choice 
of parameters A p /(2tt) = 1 GHz, Q p /{2tt) pa 10 MHz 
and fl c /(2ir) pa 100 MHz fixes the time of the pulse to 
7V pa 15 fis, much shorter than the lifetime of any Ry- 
dberg n'-S state with n' > 40. This time is also much 
shorter than the typical timescale of the atomic motion 
(milliseconds), so that we can consider the gas to be 
frozen [25l[26] for the duration of the transfer. The radius 
R' c of the sphere inside which atoms are transferred to the 
imaging state can be tuned depending on the mean impu- 
rity distance i?b and available imaging resolution. This 
can be done by selecting an appropriate pair of Rydberg 
states, to which the impurities and background atoms are 
excited. As an example, in typical experiments studying 
Rydberg-blockade effects, i?b is on the order of 10 pm [3], 
while imaging resolutions are < 3 pm. For R' c = 3 pm and 
atomic densities p ~ 10 11 cm -3 , the number of atoms in 
\B) around an impurity is N-q ~ 60. For the assumed 
imaging resolution, this corresponds to an optical depth 
~ 1 in the absorption image pixel containing the impu- 
rity, resulting in an easily observed transmission reduc- 
tion of the imaging light, caused by the single impurity. 

In conclusion, we have developed a novel method for 
spatially resolved imaging of single impurities in an ultra- 
cold gas. Our scheme provides access to direct study of 
impurity dynamics which can be straight-forwardly im- 
plemented in existing experiments. 
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During preparation of this manuscript we became 
aware of related work 27]. 
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